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Abstract 

A version of the J-matrix method for solving numerically the three-body Faddeev-Merkuriev 
differential equations is proposed. This version allows to take into account the full spectrum of the 
two-body Coulomb subsystem. As a result, a discrete analog of the Lippmann-Schwinger equation 
is obtained which allows to interpret correctly the three-body wave function in two-body domains. 
The scheme is applied to calculations of the fully resolved absolute differential cross sections for the 
He(e, 2e)He''' and He(e, 3e)He~'~'^ reactions at small energy and momentum transfers. The results 
are in good agreement with the experiment both in shape and in absolute value. 
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I. INTRODUCTION 



The few-body Coulomb problem is one of the most fundamental problems in physics and 
being yet to be exactly solved. In recent years, it has attracted numerous theoretical studies 
which involve different approaches and/or approximations. Among these is the J-matrix 
method widely used in the quantum scattering theory. It was first proposed in atomic 
physics [1] and later on (independently) in nuclear physics 3,0,0]. Within this approach 
the full Hamiltonian of the atomic system is split into two parts H = Hc + V. The operator 
He determines the discrete basis of square-integrable functions which are used for expansion 
of the many-body wave function. In atomic physics, a Laguerre basis is the most widely 
used because it provides the three-diagonal representation of a radial part of the operator 
He, and this infinite trinomial recurrence can be solved analytically. At the same time, 
the short-range interaction V is approached by its projection on a subspace of N basic 
functions. Hence, the J-matrix method supplies the exact solution of the scattering problem 
on the model potential . The J-matrix approach is proved to be an efficient and rather 
accurate numerical method. Formally and from the viewpoint of a numerical realization, 
the J-matrix method is similar to the R-matrix scattering theory (for an overview, see, for 
example, 0). It is also equivalent to the method of the so-called Coulomb-Sturm separable 
expansion (see, for instance, P] and references therein). 

The most important problem in a few-body Coulomb scattering theory is the description 
of the continuum-state wave function in terms of square-integrable functions. The method of 
pseudostates and its recent avatar, the convergent close-coupling method (CCC) replace 
the continuous energy spectrum of a selected two-body subsystem by a finite number of 
positive energy lines. However, we present here a new scheme which allows to take into 
account both the summation and integration over respectively the bound and continuum 
states of the two-body subsystem. 

The general idea guiding the J-matrix method was first formulated in jsl] for a system 
of three nuclear particles interacting with short-range potentials using an oscillator basis. 
Recently, on a base of the results of Papp et al. |10] it was extended for the case of long- 
range Coulomb potentials P]. We solve the Faddeev-Merkuriev differential equations which 
allow to formulate the boundary conditions for a component tpa of the full wave function 
\1/ = J2a in terms of the fixed number of Jacobi coordinates Xq, y^. We use the fact jlll 
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that in the two-body domain Qa, where the mutual distance Xq between particles f3 and 7 
is much smaller than the distance ja between their center-of-mass and particle a, the total 
asymptotic Hamiltonian Has looks like a direct sum of uncoupled Hamiltonians h^^ and hy^. 
As a result, the full three-body Green's function G = {E — H + iO) is connected with the 
asymptotic operator Ga = {E — h^^ — hy^ + iO) in the domain ^1^. In turn, the Green's 
function Ga is presented as an overlapping integral of the operators Gx^ and Gy^. 

Furthermore, following Ref. the component tpa is presented as a decomposition on the 
eigenfunctions of the two-body subsystem (/57) with the total charge Z^- Such a choice allows 
to factor out the long-range part of the interaction and, consequently, to describe correctly 
the asymptotic behavior of the function \1/ only in the two-body domain VL^- In this way, 
we arrive at a discrete analog of the Lippmann-Schwinger equation for the component ipa in 
VLct- This is the general scheme of the method, and the details are presented below. 

The present paper is organized as follows. A detailed theoretical formulation of the above 
general scheme is presented in Section II. In Section III, the efficiency of the presented 
numerical scheme is demonstrated in a light of comparison with experimental data. We 
calculate for a helium target the (e, 2e) triple differential cross section (TDCS) for single 
ionization accompanied with simultaneous excitation, and the (e, 3e) five-fold differential 
cross section (5DCS) for double ionization under fast electron impact. Previous calculations 
of the three-body wave function of the system (e, e,He"'"^) were carried out within the method 

n n n n 

of pseudostates 1121 113II as well as using other methods, such as the so-called 3C [U, I15| . 



approximate 6C |l4 llfil]. and C4FS |12|. Most of theoretical results exhibited a marked 



12, hj. In this work we show 



disagreement with the experimental data on absolute scale 
that the proposed method is able to describe the same experiments without any scaling 
factors. 

II. THEORY 

The Hamiltonian of a three-body system has the form 

3 

H = H, + Y.'^^^xa). (1) 

where Hq is the kinetic energy operator 

ifo = -A,„- A,„, (2) 



and 



The couple (x^, ya) stands for the set of Jacobi coordinates |ll| 
where mi are the particle masses and 



V^{x.) = ^. (3) 



nifjrfj+mjrj 



(4) 



.„-./2^^. ,.= j2,nJl . (5) 

771/3 + y V mi +1712 + 7713 J 

The interaction can be decomposed into the short- and long-range parts {Va'^'^ and 
Va \ respectively) 

with the "separation" function of the form 

Cixa,ya) = 2/ {l + exp[{xa/xoy/{l + ya/yo)]} (7) 

where u > 2. Thus, the function Va^^ decreases rather rapidly in the "true" three-body 
asymptotic domain Qq and coincides with the initial potential in the two-body asymptotic 
domain Qa {^a -C ya)- In the general case of different particles, the total wave function is 
represented as the sum of three components ipa satisfying the set of equations 

Ho + Va{xa) + V^'\x^) + V.f\x^) - ^] ^„ = -V^'^ i^p + . (8) 

For the system (e, e, He^^) particles 1 and 2 (electrons) are identical, and the solution 
reduces to the sum of two components ipi and iIj2 (see, for instance, |lSi])- They are related 
to each other as follows: 7p2 = gPu'ipi {g = +^ and g = —1 for, respectively, singlet 
and triplet spin states, P12 is the permutation operator). Taking into account the spatial 
symmetry of the total wave function, Eq. (jHJ can be reduced to a single equation [3] 

[^^0 + Viixi) + V^ixs) + V^'\x2) - E]M^i, yi) = -gV}'^ P^M^i, Yi) (9) 
for the component V'l- 
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The wave function ^ of the system (e, e,B.e'^~^) with all three particles being free can 
be represented by the following decomposition 



mo fio 

(10) 

The corresponding decomposition of the function "^^Ji^^^ of the system (e,He^) with two 
bound particles is given by 

^ L Ao /JO 

where the quantum numbers [uq io m) describe the state of the He"*" ion. 

The spatial part ipio^ component ipi in ffTUI) and can be presented in the form 

of a bispherical expansion 

#ot(x, y) = E "^ y'A"(^> ---1' y-yi' (12) 

J^if (x, y) = Y,i^m\i,\LM)YUx)Y,,{y). (13) 

m n 

I for the radial functions V'4,Ao,^a\ 
of an expansion 

^^oAo,.a(^, y) = J2 f dkai\k)^,A^)<t>tiy) (14) 
in the eigenf unctions (pk/{x) of the Hamiltonian 

/i, = -A, + — (15) 

X 

describing the subsystem (2, 3). Note that here we pave the way for nondemocratic succes- 
sive ejection of electrons from the atom. 

In (|T4j) and hereafter J dk means the summation over the discrete states and integration 
over the continuous states of the subsystem (2, 3), i.e. 

oo 

j dkai\k)ifkAx) = Yl^i^(^) + IJ dkMk, x) ai^k), (16) 

^ 

with y^^"''' and (pe{k, x) being the corresponding eigenf unctions of the Hamiltonian (fTH|l [see 

Eqs. (IA.2|1 and ()A.3j) of Appendix]. The Laguerre basis functions 0^ flj are used in (fT^ 



In turn, and in accord with j8[, we look for the radial functions ijjf ex{x, y) in the form 



with u being the basis parameter which suitable choice affects the rate of convergence of the 
numerical results. 

The functions (p[^^ and (pe{k, x) can be also decomposed on the basis functions 0^ (fTTj). 
We use the analytical expressions for the coefficients S^^j and Sne{k) of such a decomposition 
[see Eqs. ()A.5|) and ()A.6|) of Appendix]. Thus, the function V'^Aq reads 

<l = E 1^^^-^' LM) , (18) 

£, A, ra, u 

\niuX; LM) = '^M^My^M^^^ (^g) 
xy 

The coefficients Cn!f^^ are of the form 

C^(^")(E) = / dkSni{k)ai\k). (20) 



The set of algebraic equations for the coefficients al^{k) was obtained in the papers 
As a result, Eq. Q transforms to the form 



JUp)<Hk) = - E k'+p' = E. (21) 

n, Ji', 

Here J'^^^, are the elements of the J-matrix [see Eq. ()A.9|) of Appendix] corresponding to 
the operator {hy — p^) {hy = —Ay + Zu/y). The potential Zujy describes the Coulomb 
interaction of the particle 1 with the center-of-mass of the subsystem (2, 3), and values 
'^nun't' '* denote the matrix elements of the operator 



V(x, y) = Vsixs) + Vi'\x2) - ^ + gVl'^P,,, (22) 



I.e. 



Kfi'n'l'''^ = {niuX; LM\ \/(x, y) \n'i'u'X'; LM) . (23) 

Within the two-body domain where x << y — oo, the potential V^(x, y) is a short-range one. 
This circumstance allows to take into account only a finite number of terms in the r.h.s. of 
Eq. ()2H). This means that Vi^^^J ^'^ = if at least one of the indexes {n, u, n', u'} extends 
to some rather large number N. Thus, ii u > N then the r.h.s. of eq. ()2H1 turns to zero, 
and the coefficients a^^ satisfy the "free" equation 

Ju\-i{p) C-i{k) + JUp) <\k) + J^^^M a%{k) = 0. (24) 



Now we use the expression 



St^^\p) = --<Su< \{p) Cl^'xip), ^< = min {z/, z/'} , z/> = max {u, u'} 



P 



(25) 



for the matrix elements [^^'•^-'(p)] of the radial Green's function. The functions C^^{p) are 
determined by Eq. ()A.14|) of Appendix. Since the infinite matrix [^'"'^^•'(p)] formally is the 
inverse matrix to the infinite J-matrix [^/'^(p)] [2^, i. 



I.e. 



-[g''^\p)] [j\p)] = i, 

the set of equations (PT|) can be rewritten in the following form 

ai^{k) = S^XoiPo)S{e\){eoXo)S{k - ko) + 



(26) 



(27) 



+ yuu' {P)'^n'£[k) V„,^,^„^„ G„„^„ {h}, Z/-U, 1,.... 

n',u',n",u" 

e",x" 

Inserting (j27|) into (j2(Jj) and taking into account the boundary conditions for the functions 
(fTUj) and (EH)' obtain the set of equations for the expansion coefficients 

Cni^^^ in (dHl). In the case, where all particles are asymptotically free, it takes the form 



Cnu^\E) = X) (eo \o) ^nioi^o) ^uXo{Po) + 

+ E Ijdk S^eik) S^,,{k) g'J-^ (p)] 'i' CtZ''""^ (E) 



n',u',n",u" 
i",X" 



(28) 



In the case of two asymptotically bound particles the set of equations takes the form 



a^(")(E) = 6^,x)ieoXo)SZ'SuxM+ 

+ E \jdkSr.eik)Sr,eik)g;:i-\p) 



-(no) 



n',u',n" v" 
l",X" 



..L(£A)(£" X") ^L(l" A")/px 
^n'v'n"v" ^n"u" \^ I 



(29) 



The set of equations and (PU)) is in fact the discrete analog of an integral equation of 
the Lippmann-Schwinger type, which was obtained in [IQ,]. It is easy to see that the integral 
in square brackets of ()28|) coincides with the matrix elements G^^l j^, of the asymptotic 
three-body Green's function G^^"^ in the two-body domain. It has the form of a contour 



overlap integral of the two-body Green's functions Gx and Gy ^ 

G^-\E) = ^JdS Gi-) (^)G(-) (E-S). 



(30) 



Here the contour C surrounds the spectrum of the operator hx in a anticlockwise direction. 

The integrals in the r.h.s. of Eq. ()28|) can be calculated directly. Note that the integrand 
has the following poles at the points kj = E + r]j {rjj = — Zii/[2(A + j + 1)]) which 
correspond to the discrete spectrum of the Hamiltonian hy and have a densening point at 
fcoo = V^- The presence of the matrix element of the radial Green's function Q^\^, "* (p) in 
the integrand means that the poles have to be rounded in a clockwise direction (Fig. 1). 

Let us note that the poles can be allocated more uniformly along the contour of in- 
tegration if one carries out the transformation of variable. For example, if one puts 
q = —Zu/C^Vk^ — E) then the poles are disposed at the points = A + 1 + j. Thus, 
the integral under consideration can be performed as an infinite sum of integrals Ij along 
segments Lj, in which the poles are localized, and integral le with a smooth integrand, 
which is calculated along the segment Lg disposing to the right of the segment Lq. Using 
the Sokhotsky formula we obtain 

I, = ^idkSni{k)Sn'e{k)gt^^\p)-2iSne{k,)Sn'e{kj)Resg^i-\p) = 

Lj 

= I idkSneik) Sr,,{k) g^i-\p) + iSUkj)Sn'e{k,) S^^, 

For E < the matrix elements G^,^).l^, can be calculated without difficulties because 
the number of poles is finite and (usually) not numerous. For > we must calculate an 
infinite number of integrals in (j31|) . In this particular case we calculate the overlap integral 
(|3(J|) being written now for the matrix elements 

gI'XUe) = ^J d£gl^-\Xs^)gti-\X8;^l (32) 
c 

with £q = ^E, and use the method of Shakeshaft The poles and cuts of the integrand 
in (jH^ are shown in Fig. 2 (here e is not infinitesimal for visualization) . To simplify 
calculations we rotate the contour C by a positive angle (p, and the new contour C = Soe^'^t 
is depicted by the dotted line. 

The integral (j32j) is calculated numerically along the contour C. To explain details of 
integration we subdivide the path of integration conventionally into two parts: C[ {t < 0) 
and C2 (t > 0). We also introduce for brevity two variables: p = y/So — S and k = v^J^J~+T. 
Note that arg(p) > for t < 0, whereas Cl^{p) and g^\^,\p) are defined at — vr < arg(p) < 0. 
To define g^^ip) on C[ we use the analytic continuation [19^ Cl^{p) = Cl^x{p) — 2iSyx{p)^ 
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(31 



or 



,A(+), 



1i 



(33) 



At the same time, the contour C[ passes in the domain of definition of the function Q^^^){}i). 

Reasoning by analogy allows us to integrate along the contour C^- Finally we obtain in 
the limit e — > 





-oo 
oo 



+ 



+ / OH^) + TSniik) Sr^nik) gyip) 



Calculating the second term in (jHS)) we accept in accordance with that 

\r{£ + l + it)f ^T{e + l + it)T{i + 1 - it), 



(34) 



(35) 



as well as calculating ^ we choose the minimal meaning of the argument ^ at the range 
—71 < arg(^) < 71. 

III. RESULTS AND DISCUSSION 



To illustrate the efficiency of the presented numerical scheme, calculations of the triple 
differential cross section (TDCS) for (e, 2e) single ionization and five-fold differential cross 
section (5DCS) for (e, 3e) double-ionization reactions on the helium atom in a singlet state 
were performed. If a fast projectile electron of energy about several kiloelectronvolts transfers 
to the atom relatively small amounts of energy and momentum, the respective four-body 
problem can be sufficiently simplified by examining only the first Born approximation in the 
interaction of the projectile electron with the atom. The calculations were performed in the 
limit ms =— > oo; that is, Xi = \/2t2 and yi = -\/2i"i (atomic units are used me = e = h = 1). 

The triple differential cross section (TDCS) of He(e, 2e)He"'" reaction when the residual 
He^ ion remains in an exited state is written in the form 



a, 



(3) 
no 



'no 



dVtsdEi dVti 



to m 



(-) 



exp(zQri) + exp(?Qr2) - 2 \^q) 



(36) 



where (i?j,pj), (i?s,Ps), and (-E'i,ki) are the energy and momentum of, respectively, the 
incident (fast), the scattered (fast), and the ejected (slow) electron; Q = — is the 
transferred momentum. 
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The five-fold differential cross section (5DCS) of He(e, 3e)He reaction is given by 

- ,n.,EXm<ia. - W exp(.Qr,) + exp(.Qr,) - 2 ^ • 

(37) 

where ki) and are the energies and momenta of both ejected slow electrons. 

The final-state wave functions (HH and ^^j^^ (HH) are obtained using the method 
described in the previous Section. We have restricted ourselves to the maximum value of the 
total orbital angular momentum, L^ax = 2, and £, A < 3. The number of basis functions 
fll7|) for each Jacobi coordinate x and y is taken as = 21. The basis parameter u fll7|) 
affects the rate of convergence of the numerical results, and its optimum value turned out 
to be M = 0.3. We also use the following parameters of the separation function v = 2.5; 
xq = 0.9; yo = 6. 

The helium ground-state wave function \E'o is obtained as a result of diagonalization of 
the matrix which was calculated in the basis (fT^ . The orthonormal functions 



2u 



{2urY+'e-'''Li'+\2ur) (38) 



are used here instead of functions (fTTj). Here we put imax = 3 and rimax = ^max = 15. 
Choosing the basis parameter uq = 1.193 yields the value Eq = —2.903256 for the ground- 
state energy. 

The angular distributions of the slow electron in the case of (e, 2e) ionization-excitation 
reaction are presented in Figs. 3-5, 6i being its in-plane emission angle with respect to 
the vector pj. The displayed experimental data, which correspond to the general kinematic 
conditions uq = 2, Es = 5500 eV, and three particular cases Ei = 5 eV and 9s = 0.35° (Fig. 
3), Ei = 10 eV and 9^ = 0.32° (Fig.4l^ind Ei = 75 eV and 9s = 1° (Fig. 5), where 9s is the 
scattering angle, are borrowed from (22j |. 

The results for a^^^ are presented in Figs. 6 and 7. The in-plane angle 9i of one of the 
two slow electrons is fixed, while the in-plane angle 92 of the other slow electron varies. The 
energy of the scattered electron Es = 5500 eV and its in-plane angle 9s = 0.45° is also fixed 
in all experiments. The energies of the slow electrons are Ei = E2 = 10 eV (Fig. 6) and 
El = E2 = 4: eV (Fig. 7). One can see that our results quite satisfactorily agree with the 
experimental distributions both in shape and in absolute value. This agreement in absolute 
scale favorably distinguishes our calculations from that obtained earlier in by the 
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method of pseudostates. Considerable scaling factors were needed there to compare theory 
and experiment. 

Let us discuss this success of our treatment in more detail. The exact final state wave 
ftmction |^(^^(ki,k2) > must be normalized 

< ¥-\K, k',)\¥-\k^, k2) >= 5(ki - k[)5{k2 - K) + 5(ki - k^)(5(k2 - k;) 

The symmetrized sum of plane waves obviously possesses this property, as well as do 



1$ 



Coulomb waves. Recent calculations with a 3C function 

vl>(-)(ki,k2;ri,r2) = ^^^^ e-^'^^-^^-v^{k„r,)^^{k,,r,)if^,{ku,r,2) (39) 

and some variational ground-state helium functions have demonstrated striking agreement 
with the experiment, although the published evidences that this function is normalized are 
unknown to us. 

It is also obvious that the function '^^^^ being decomposed into a finite sum of square- 
integrable functions can never be normalized to 6 function. Our numerical scheme allows 
to conclude that the asymptotic bound of the function \E''-~'' is a product of two Coulomb 
functions in the discrete representation. Presumably, it is this property which provides 
much better normalization conditions for the continuum wave function than in the case of 
the pseudostates' approach. 

As a summary, we formulate two main conclusions 

• The proposed numerical scheme and calculations demonstrate the importance of ac- 
counting for the whole two-body continuum spectrum. The method of pseudostates 
which replaces the continuum by a finite number of states with positive energies, 
seems to face with serious difficulties as to the magnitude of the calculated differential 
cross sections, especially when the resultiiig final-state wave function is applied to the 
calculation of (e, 3e) matrix elements |l^ Ifsj]. 

• We observe a lack of coincidence between the theory and experiment in some kine- 
matical cases. This, perhaps, testifies to a necessity of taking into account the correct 
behavior of the function \1''^^) in a "true" three-body asymptotic region. The simple 
3C model (jSHl) partially accounts for such a behavior. 

Finally we can conclude that the presented method based on the J-matrix approach 
allows to formulate the effective numerical scheme for applications in atomic physics. 
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APPENDIX 



The Hamiltonian 



H, = -^^+l{l + l)/r-' + 2Z/r 



has the following wave functions of the bound spectrum {Z < 0) 



l3) 



(r) = (2Kj-r)^+ie-''^" 



\/-Z (j + 1)2^+1 



(j+£ + l)(2£+l) 



(A.l) 



(A.2) 



and continuous spectrum 



ipi{k, r) = -{2kr) 



2^--' (2£ + l)! 

They can be expanded in the Laguerre basis functions [l| 

0lM = [{n + l\,,^,)]-''\2urY^'e-^^Ll'^\2ur) 

and the expansion coefficients S'^l and Sniik) are given by |l, 19| 

Si^J = (-1) 



iFi{i+l + it; 2i + 2; -2ikr). (A.3) 



(A.4) 



^+1 



(i+^+l)(2^+l)! 



X 



(A.5) 



and 



Sniik) = \ ^{n + 1)2^+1 (2 sin 0^+^ e^-^/^ C 



|r(£ + l + zt)| 



(A.6) 



(2£ + l)! " 
X (-0" 2F1 (-n, £ + 1 + it; 2£ + 2; 1 - T') • 
Here = -Z/U + i + l),t = Z/k, and ^ = e^^ = {tu - k)/{iu + A;). 

It is proved 1| that the functions S^^} , Sn.e are the regular solutions of the infinite trinomial 
recurrence which is a discrete analog of the Schrodinger equation 

Jn,n-lik)d^-l + Jn,nik)dn + Jln+^drr+l = (n = 1, 2, . . .) (A.7) 

with the initial condition 

jloik)do + Jli{k)di = 0. (A.8) 

In flA.7|l J'^ „/(fc) are the elements of three-diagonal matrix of the operator {Hi — k^) (i.e. 
J-matrix) calculated with the basis functions 0^ 

jUk) = ^{n + i+l) + 2Z, 

JL-Ak) = ^v/^^TTWTT), (A.9) 



JLAk) = ^V(n + l)(n + 2£ + 2). 
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The coefficient functions S^j ^ Sni satisfy the normahzation conditions 



n. n'=0 



(A.IO) 



n, n'=0 



where Ql^^i are the elements of three-diagonal matrix and are given by the overlapping 
integrals 



Qnn' 



(A.ll) 



QL = iin + i+l). 
The corresponding completeness condition takes the form 



(A.12) 



j=0 



The irregular solution of Eq. ()A.7|) can be written in the form jl. Il9| 



-0 



±(n+l) 



2F1 (-£ ± it, n + 1; n + £ + 2 ± 2t; ^ 



±2\ 



(A.13) 



(A. 14) 



r(n + £ + 2 ± it) 

The function Cl^g^{k) (^Cl^^ik)^ is determined at Im(fc) > [Im(A;) < 0] of the compl ex p lane 
k. The analytical continuation can be done with the help of the following relation |l9|: 



ci;',\k) = c'-,>{k) + 2tSni{k). 



(A.15) 
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FIG. 1: The integration contour of the integral in parentheses in Eq. ()28() . The poles of the 
integrand are shown as closed circles. 
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FIG. 2: The path of integration of the convolution integral l|H2|) . The poles which are associated 
with the discrete spectrum of the Hamiltonian hx are depicted as closed circles. —Sq + ie denotes 
the initial point of the corresponding unitary branch cut. Open circles and — ie are the same 
for the Hamiltonian hy. 
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FIG. 3: TDCS for the ionization-excitation reaction He(e, 2e)iie~^ when the hehum ion is left in 
the state no = 2. The fast scattered electron energy is Eg = 5500 eV, while Ei = 5 eV for the 
slow ejected electron. The scattering angle of the fast incident electron is fixed Og = 0.35", and 
the angle of the ejected electron 9i varies relative to the incident electron direction. All electron 
velocities are disposed on the same plane. The experimental data are borrowed from 22 1. 
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FIG. 4: The same as in Fig. 3, except Ei = 10 eV and Os = 0.32°. 
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FIG. 5: The same as in Fig. 3, except Ei = 75 eV and 6s = 1°. 
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FIG. 6: The fully resolved five-fold differential cross section 5DCS of the electron-impact double 
ionization reaction He(e, 3e)He^'''. The energy of the scattered electron is Eg = 5500 eV, and the 
energies of the slow ejected electrons are Ei = E2 = 10 eV. The scattering angle of the fast incident 
electron is fixed Og = 0.45°, and the angles of the ejected electrons are 61 and 62^ where one angle is 
fixed while the other angle varies relative to the incident electron direction. All electron velocities 

n 

are disposed on the same plane. The absolute measurements are borrowed from |17l |. 
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FIG. 7: The same as in Fig. 6, except the ejection energies are lowered to Ei = E2 = 4: eV. 
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